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Abstract 

We compare some recent computations of the entanglement of forma- 
tion in quantum information theory and of the entropy of a subalgebra in 
quantum ergodic theory. Both notions require optimization over decom- 
positions of quantum states. We show that both functionals are strongly 
related for some highly symmetric density matrices. Indeed, for certain 
interesting regions the entanglement of formation can be expressed by 
the entropy of a commuting subalgebra, and the corresponding optimal 
decompositions can be obtained one from the other. We discuss the pres- 
ence of broken symmetries in relation with the structure of the optimal 
decompositions. 

1. INTRODUCTION 

Entanglement, always one of the most intriguing among quantum marvels, 
has lately become a powerful resource in prospective quantum information tech- 
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nologies |T|; measuring the entanglement content of states of multipartite quan- 
tum systems is thus of great practical importance. If a bipartite system A + B is 
described by a density matrix pab, the so-called entanglement of formation || 
is measured by 

E(pab) ■= inf{Xl ^jSiTiB-Kj) : Pab = X] V 1 ".?} • (*) 

j j 

In the above expression, S(p) := — Trplogp denotes the von Neumann entropy 
of the state obtained by partial trace over B and the infimum is computed over 
all possible decompositions of p as convexly linear combinations, that is Aj > 0, 
Aj = 1, of one-dimensional projections ttj of A + B. In the following we call 
such decompositions extremal convex decompositions of p to be distinguished 
from generic convex decompositions into mixed states. 

When pab = \^?ab)($?ab\, the entanglement of formation gives the asymp- 
totic ratio between the number of singlet states necessary to construct N ^ 1 
copies of pab 0- In the following, we will compare the entanglement of forma- 
tion with a particular case of a more general quantity, the "entanglement with 
respect to a subalgebra" or "entanglement", for short. This latter concept is 
related to the so-called "entropy of a subalgebra" A contained in a reference 
algebra A4, relative to a state p on M. ||, 

H p {A) := S(p \A) - inf X 3 S( Pj \A) : p = £ X jPj } . (2) 

j j 

In the above expression, the infimum is calculated over all convexly linear de- 
compositions of p into other states on M.. It plays a key role in extending the 
classical dynamical entropy of Kolmogorov to quantum systems |5], || [7]]. The 
entanglement of formation (|I|) can be considered a special case of (0). 

We shall call "optimal" those decompositions achieving the extremum in ([I]) 
and (0). Calculating either E(pab) or H P (A) is particularly complicated. The 
problem has been completely solved for the entanglement of formation if H,a = 
H B = C 2 [|, and for the entropy of a subalgebra if M = M 2 (C) 0, |, [K|. 
So far, all other available results concern states pab and p that are highly 
symmetric, isotropic in [ H|] , respectively permutation- invariant in |TB| . 

In this paper we will discuss the previously mentioned results by comparing 
the two notions of entanglement sketched above. We show, that some of them 
are one-to-one related. To do so, we shall focus on the structure of optimal 
decompositions in relation to the symmetries existing in the problem and show 
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possible ways of breaking them. These symmetries form a group G and leave 
invariant both the state p and, as a set, the subalgebra A. Given extremal 
optimal decompositions, the G-orbits of each of their pure states consist of 
optimal decomposers, too. We will study the dependence of either entanglements 
upon the number of different orbits. 

2. ENTANGLEMENT 

In the following, we shall consider quantum systems described by algebras 
of operators, A4, acting on finite or infinite dimensional Hilbert spaces 7i, with 
states, M. 3 X i — > Tr(pX), represented by density matrices which we shall 
denote by greek letters. 

Definition 2.1 Given a finite dimensional subalgebra A C Ai, we define the 
entanglement of the state p with respect to A by 

E(p; M,A):= inf {£ \ 3 S{p 3 \A): p = £ \ jPj ] , (3) 

3 3 

where p = J2j ^jPj runs through all convexly linear decompositions of p with 
states of M., and S( P j \A) is the von Neumann entropy of the state pj restricted 
to the subalgebra A. The entanglement (^) is convex as a function of p. 

Remarks 2.1 

(i) The entanglement (j^) is a convex functional over the states: 

E(J2 Wi\ M,A)<J2 mE( P j\ M,A) , £ Hj = 1 , Hj > . (4) 

3 3 3 

This follows by choosing optimal decompositions for the pj's, which together 
provide a decomposition, not necessarily optimal, for J2j PjPj- 

(ii) The entanglement of formation (Q) is the entanglement (|3]) with A, respec- 
tively B, the algebra of observables of the system A, respectively B, M. = A®B 
and p A B \A = Tt b Pab- 

(in) The entanglement @ is related with the entropy of a subalgebra (||) by 
E(pab) = S (pab \A ® 1 B ) - H PAB (A ® 1 B ) • (5) 

Indeed, as we shall see below in Proposition 2.1, the infimum is achieved at de- 
compositions using pure states of M. only, and it enjoys some further remarkable 
properties. 

The quantity in ([5]) and some techniques 0, 14 1 that were developed for 
computing (H), have recently been used to attack the question whether the 
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entanglement of formation is additive |15|]. Among them, a useful result is con- 



tained in the following proposition. The idea is in [y| and, slightly extended, 
in [[T^]. We include a proof for the benefit of the reader. 

Proposition 2.1 If the algebra M. is finite dimensional then 



• the entanglement Eyp;A4,Aj is achieved at certain extremal convex de- 
compositions p = J2j ^j^ji <\j > which saturate (|3|). Such decompositions 
are called optimal. Every pure state, ir, which appears in at least one op- 
timal decomposition of p is called p-optimal or an optimal decomposers of 
P- 

• For every p there is an optimal decomposition with a length not exceeding 
the linear dimension of A4. 

• The functional . ; At, Aj is convexly linear on the convex hull TZ{p) of 
all p-optimal pure states: Let be u = J^iC^i^i: a i > 0, J2i a i — 1 an y ex- 
tremal convex decomposition where the iTj are some optimal decomposers 
of p. Then 

E(u;M,A)=^2^S(n l \A) . (6) 

i 

Proof: Any mixed state p can be decomposed and, since the von Neumann 
entropy is concave on convex combinations, mixed states cannot improve (§) 
with respect to pure states. If Ai is d dimensional, compactness of the state 
space, extremality and compactness of the set of pure states ensure by a the- 
orem of Caratheodory that we need not less than d and not more than d 2 
decomposers JT0], |TjJ. Because of convexity (|]), the functional E(^.;Ai,Aj is 
the supremum over affine functionals. Thus, for every p there are functionals I 
such that E^p;Ai,Aj = i(p), while, on generic states a, E(a;Ai,Aj > £(a). 
Given an optimal decomposition p = J2j ^j^j h follows 



e(j>-,m,a) = Y, x i E fa> M > A ) 

j 

> J2 X A^)=e(p)=E(p;M,A) . (7) 

i 

Since equality must hold in (|7|) and because Aj > 0, while eOkj; Ai, A) > £(p) 
by assumption, we conclude E^jif, Ai, A^j = ((ttj) for all j. With u e 7£(p), 
let us now fix this affine functional i and consider the extremal decomposition 
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uj = &kir' k such that all the 7r| are optimal decomposers of p. By convexity 
and the preceding argument we deduce 

k k 

However, i{uo) < e(u; Ai, Aj by our choice of £, and equality holds in (||). 
Thus, E(^-; Ai, Aj is convexly linear on TZ(p). ■ 

Definition 2.2 We shall call the convex hull TZ(p) of the optimal decomposers 
of p a leaf with respect to the entanglement E(p; Ai,Aj. Then, the state space 
appears as covered by leaves, and the entanglement itself is convexly linear 
above every leaf. That effect is referred to as the roof property of Ey ■ ; Ai, Aj, 
|10fl , i.e. E(^ ■ ; Ai,Aj is a convex roof. 

Definition 2.3 Given p on Ai, we shall call a group G a symmetry group with 
respect to E(p; Ai, A), if for all g G G there exists a linear map 7 9 : Ai i — > .M 
such that the state and the subalgebra ^4 (as a set) are left invariant by j g , 
Namely, l* g [p] = p, where J g [p](m) = Tr(pj g (m)). 

Proposition 2.2 If G is a symmetry group with respect to E(p; Ai, A), the 
leaf 7Z{p) is G-invariant as a set. In particular, the action of G permutes the 
optimal decomposers of p. 

Proof: Let p = ^jPj be an optimal decomposition with respect to 

E(p; Ai, A). Then, since 7^[p] = p and j(A) = A for g G G, the decompo- 
sition p = ^jlg(Pj) is also optimal. Therefore, its leaf TZ(p) must contain 
both the p/s and the 7g(pj)'s. ■ 

Based on the previous two propositions, the entropy H p (A) has explicitly 
been computed in the following cases, 

Case 1. B HI Let M be the ful1 2x2 matrix algebra M 2 (C), .A the 



subalgebra of all 2 x 2 matrices diagonal with respect to a given basis |1), |2), 
/ a b \ 

and p = [ }, < a < 1, \b\ 2 < a(l — a), any density matrix. 



Case 2. [13| Let .M = M 3 (C), A the subalgebra of all 3 x 3 diagonal matrices 
with respect to the basis |1), |2), |3) and 

/ 1 x x\ 

, -1/2 < x < 1 , (9) 



x 1 a; 
\i i 1 J 
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any density matrix invariant under the group of permutations of {1, 2, 3}. 

For future comparison with the entanglement of formation of isotropic states 
of d- dimensional bipartite systems studied in [11 1, we fix an ortho normal basis 
\j) G C d and consider the group G of permutations of {1, 2, . . . , d}. It turns out 
that any G-invariant density matrix p(x) over M. = M^(C) can be written as 

1 — F 



Pf 



d- 1 



(10) 



where = —j= ^ |j) and F is the fidelity parameter 
V d j=1 

^ /it / \ \ i \ (d — l)x + 1 
< F := (iP\p(x)\il;) = K ^ < 1 . 

Setting s(t) := — tlogt, we have, 

Case 1. For all p, the optimal decompositions are 

p = X\wi)(w 1 \ + (1 - X)\w 2 )(w 2 \ 



z l\ 



Zl 

Z 2. 



\w 2 ) 



ZiZ 2 . 



1-4|6| 2 )/2 = 1-|^ 2 | 2 , A 



1 + 



2a- 1 
1 -4|6| 2 , 



(12) 
(13) 

.(14) 



The corresponding entanglement is E (p; M 2 (C) , A) = s( | | 2 ) + s( 1 2^2 1 2 ) • 

If p = pf is permutation-invariant, that is, if a = 1/2, b = x/2 F — (1 + x)/2, 

the entanglement reads 



E(p F -M 2 (C),A) = s 



1 + 2JF(1-F) 



+ s 



1-2JF(1-F) 



(15) 



Case 2. Given the group G of permutations of {1,2,3}, let V, V 2 implement 
unitarily the subgroup Gq of cyclic permutations. Then, any G-invariant state 
Pf can be written 



p F = -\w)(w\ + -V^iwlV- 1 + -V 2 \w)(w\V- 2 
3 3 3 



(16) 



where 



\w) 



I a + 2b cos 8 \ 

a - 2bcos(6 - tt/3) 
K a - 26cos(6> + tt/3) J 



X 3F 



;i - f) . (17) 
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The structure of optimal decompositions depends on the convexity of 



S(F) := mmV S (h(F; 



6»e[0,27r 



For F > F* := (2x* + l)/3, x* = —0.4150234, the minimum is achieved at a 
single extremal Go-orbit generated by the vectors 

/ y/F- 



\w) = - 
1 1 3 



( a + 2b\ 

a — b 
\ a-b j 



1 

7! 



2(1 



F 
F 



(l-F)/2 



(19) 



-F)/2J 

For each < F < F*, there are two different orbit-generating vectors, \w±(F)), 
whose Go-orbits provide different optimal decomposers for flT8|), and which form 
together one orbit of the full permutation group G. They are 



\w ± (F)) 



( 



\ 



a 
\ a 



(20) 



a + 2b cos ap 
26cos(7r/3 =f ol f 
26cos(7r/3 ± ap) ) 

where the angle ap varies with < F < F*. 

Finally, for F = 0, ap = —tt/6, the minimum is achieved again at a single 

G-orbit containing the vector, |w ) = —^=(1,0,-1). As the 6 vectors coincide 

v2 

pairwise up to a sign, the states form a single optimal decomposition of length 
3. 

In HOI), it is shown that the above vectors give optimal decompositions as 
long the function S(F) is convex. Numerically, this is the case for all F < 8/9. 
The corresponding entanglement is 



E(p F ;M 3 (C),A) 



2 -F + 2J2F(1 - F) 



+ 



2s(- 



1 + F — 2y2F(l — F). 

6 ■ 



(21) 



for fidelities F* < F < 8/9. For F = the entanglement equals log 2. We have 
only numerical results within the interval < F < F*, [|TJ]], reflecting that the 
exact dependence of the angle ap in ( |20"D as a function of F is unknown. 
Remark 2.2 Permutation-invariant states as in (pTOj ) can be written as aver- 
ages over the unitaries U„ implementing the permutation group G, 

1 EC^WIEk, (22) 



PF 
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if and only if |(^>|</>)| 2 = F, where \ip) is the vector in flTTI). Necessity comes 
from the fact that U n \ip) = Sufficiency: The identity 1 and form a 

basis for all possible contributions to the averages (^). 

In view of the structure of the optimal decomposers discussed above, we 
introduce a notion of regularity with respect to a subgroup of a symmetry 
group, as follows. 

Definition 2.4 Given a symmetry group G with respect to E(p; Ai, A), we 
shall call a leaf TZ(p) regular of order n with respect to a subgroup H C G, if 
there exist n pure states pj G 7Z(p) such that ■Jhlf'j] = Pj f° r a ^ h E H, whereas 
the convex span of the orbits |7g[Pj]} eG is the whole of TZ(p). 

We illustrate the previous definitions with some examples. 

Example 2.1 Let M. be a full d x d matrix algebra on C d and A C M. di- 
agonal with respect to a chosen orthonormal basis {\j)} d = i in C d . Let p be a 
symmetric density matrix, (j'|p|/c) = (k\p\j). Then, with respect to the chosen 
representation, the transposition T respects both the state and the subalgebra 
A. Also, 1Z(p) is regular with respect to G = H = {id, T}, the order of regu- 
larity depending on the state p. In fact, let n = \4>)(ip\ G TZ(p), then, because 
of Proposition 2.2, T(n) = ir' = G 7£(p), too. If n ^ n', we may con- 

sider the state u = tt/2 + tt'/2. which, by Proposition 2.1, is already optimally 
decomposed. Also, 

E(u; M,A) = S(n \A) = S(u \A) . (23) 

Instead, the decomposition 

u = uri-r n n+ _j vvrir / / n _ ^ w \y eTe (24) 

need not be optimal. However, the concavity of the von Neumann entropy yields 

S(w;.M,.A) < S{n + \A) 

+ l-^((^)) g(7r _ M) , (26) 



It thus follows from ( pB"D that 7r = 7r± ["^4, whence the components ip(i), ip'(i) 
of t/) and ip' must coincide apart from an overall phase. Thus, n = tt' and the 
T-symmetry cannot be broken. 
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Example 2.2 Let M. = A® B, with A and B isomorphic and a : A i — ► $ 
the algebraic exchange of the two of them. If p is a state on M. such that 
p o (er -1 (g) cr) = p, in general, a -1 (g> cr does not belong to any subgroup of 
regularity of p; indeed, if A (and thus B) is a <i-dimensional matrix algebra and 
{\£)} is an orthonormal basis in the corresponding Hilbert space TCa (and thus 
also in Hb), the density matrix 

p^:=i|l><l|®|2><2| + ~|2)<2|®|1)<1| , (27) 

is such that Tr (p^- 1 <g> cr)(X <g> F)) = Tr(p(X ® F)). Also, p^B is already 
optimally decomposed, E(pab', A, Ai) = is achieved with the decomposers 
|l)(l|(g>|2)(2| and |2)(2|<8>|1)(1|, which, however, are not invariant under <j^ 1 ®)a. 



Example 2.3 Let M. = A ® B, with A and B both d x d full matrix alge- 
bras. We fix the same orthonormal basis in both Hilbert spaces 1~La,b and 
consider the one-parameter group U of unitaries 

U t :=^^- h ^\j)(j\ ®\k)(k\ . (28) 

The density matrix p AB := J2j,kRjk \j)(k\ ® \j){k\, , -R = [R jk \ > 0, Tri? = 1, is 
[/-invariant; moreover, ^/pAB = J2j,k(yR)jk\ \j){k\®\j){k\, so that the operators 
y/PABM^/p AB , M E /A, have the same matrix structure as pab- Choosing 
positive Mj > 0, j E J, such that Mj = 1, pab decomposes into 

Since it is also true that every mixed state p on .M can be written as in (|29| ) by 
means of a suitable positive Mj, ( p9|) indeed exhausts all possible decompositions 
of p^B. Thus, the decomposers 7Tj of pab which are optimal with respect to 
E(pab] M.,A), have the same structure of pab and are then [/-invariant. Hence, 
the group U is a group of symmetries of pab with respect to entanglement and 
the leaf IZ(pab) is regular with respect to H = U, its order depending on which 
further symmetries are enjoyed by pab- 

Example 2.4 Let M. = M 2 (C),A as in Case 1, and pp a permutation- 
invariant state. The leaf TZ(pp) is the orbit of the group G of permutations 
of {1,2}. This follows from the form of the optimal vectors ( |P2"D in such a case: 
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\w x ) = y z J, \w 2 ) = ^ , with z lj2 = y 1/2(1 ± 2yF(l — F)). It is regular 

of order 1 with respect to rotations with elements from A. 

Example 2.5 Let M. = M d (C) and pp a permutation-invariant state. Then, 
for F* < F and F belonging to the convexity region of S(F) in flI8|), the 
structure of the optimal vectors (|l^) ensures that the leaf TZ(pf) is regular 
of order 1 for the subgroup H of permutations {2,3} i— > {3,2}. However, at 
the point F = F* such a ii-invariant vector bifurcates into the two optimal 
ones (|20|). Thus regularity with respect to the subgroup H is broken and remains 
broken for < F < F*. At F = optimal vector states of different Go orbits 
degenerate pairwise into a single one, and one of them is if- invariant, while the 
corresponding vector changes its sign. 

In the last two examples, for all F when d = 2, and for F greater than 
the bifurcation values F* in the convexity region of S(F) in (|T3|), when d — 3, 
the leaf IZ(pp) of a permutation-invariant pp is generated by the orbit under 
the subgroup G of cyclic permutations V^\w), j = 0,1,2. The vector \w) is 
invariant under a unique transposition out of G. This structure is indeed more 
general as will be showed in the next two propositions. 

Proposition 2.3 Let A C M. = M d (C) be chosen as in Example 2.1 and the 
density matrix pp be invariant with respect to the permutation group G. If the 
leaf TZ(pp) with respect to A is generated by exactly one Go-orbit of a normal- 
ized vector state \w) G C d , with G C G the subgroup of cyclic permutations, 
then the entanglement is 

E(p F ;M d (C,A) = s{pp) + (d-l)s( ] j^-) (30) 



Pf 



y/F +v /(rf-l)(l-F) 
d 



(31) 



Remarks 2.3 

(i) The assumption of the previous proposition amounts to ask TZ(pp) to be 
regular of order 1 with respect to the subgroup H C G of permutations on 
{2, 3, ■ • • , d}. Indeed, the leaf is G-invariant, so that the d states \<pj) = V^\w), 
j = 0, 1, . . . , d — 1, obtained via cyclic permutations, must be invariant under 
the remaining (d — 1)! permutations This is possible only if d — 1 of the d 
components of the optimal vector \w) are equal. 
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(ii) If \w) has three different components, then the decompositions ( p^) con- 
tain at least d(d — 1) different terms. 

(Hi) In section 3 we will show that, upon identification of pp with the quantity 
7(F) in [11 , the entanglement of formation calculated there is given by (|31|) and 



(|30|) in a range F** > F > 1/d. The upper limit F** is a particular bifurcation 
point which was discovered in JTT] and that will be reinterpreted accordingly 
within the framework of this work. 

^ d-l 

Proof: By hypothesis, pp = — ^ V 3 \w)(w\V~ :! is an optimal decomposition 

d j=o 

with entanglement 

E(p F ;M d (C),A) = Y,s (\(j\w)\ 2 ) . (32) 
i=i 

Also, taking into account Remark 2.2 and 2.3, and decomposing 

d 

\w) = VF|V> + eVT^FlO = a\l) + , 
where £ is a pure phase, it follows that \w±) = (Vd\l) — \ip))/\fd — 1 and 

H = -j= [(VF + ey/(l-F)(d-l))\l) + (VF- ef^) E \j) 
With ^ := 2TZe(e), the right hand side of ( p2|) reads 



i=2 



d- 1 



5(0 = s(p(0) + (d-l)a( 

F + (1 - F)(d - 1) + e V /i? ( 1 - F)(d - 1) 

p(0 = 5 • 

It achieves its minimum at the maximum value of p that is for e = 1, from 
which the result follows. Indeed, as we show below, \w) must be real. If remark 
2.3(i) applies we always get a local extremum. Either by direct calculation or 



relying on |T3[ one concludes e = 1. 



We now relax the hypothesis of the previous proposition and allow for more 
than one Go-orbit to be optimal for the entanglement of pp with respect to the 
subalgebra A, that is we allow the leaf IZ(pp) to be generated by more than 
one Go-orbit. 

Proposition 2.4 Let A C Ai = M d [C) be chosen as in Example 2.1. If the 
density matrix pp is invariant with respect to the permutation group G and its 
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\w) = —j= \k) in which case F = 1 and pf = \ip){ip\] 



entanglement with respect to A can be achieved at an optimal decompositions 
consisting of one Go-orbits of normalized vector states \w) G C d , with Go C G 
the subgroup of cyclic permutations, then we have three possibilities 

d 

• \w) is real with 1 component equal to ai and d — 1 real components all 
equal to ^ a\\ 

• \w) is real with 2 components a\ ^ 03 and d — 2 components all equal to 
03 different from both a\ and (Z2. 

To prove the result we need a preliminary 

Lemma 2.1 The vector \w) whose Go-orbit is optimal can be chosen real. 

Proof: Let v^, k = 1, 2, . . . , d, be the components of \w) with respect to the 

1 d 

chosen orthonormal basis {\k)} and = —j= ^ The assumption is that 

■y d-1 

Pf = — ^ V"'|w)(w|V r_J ; from normalization it follows that the components of 

« j=0 
\w) must satisfy 

d d 2 d 

^2 \wk\ 2 = 1 , I XI w fc| = 1 _ w *£ w i = dF ■ ( 33 ) 

k=l k=l t^k=l 

Further, in order to implement optimality and achieve E(pp] Ai, A), we mini- 
mize 

d d 
S(w, A,/i) := - X |wfc| 2 log \w k \ 2 + A X \w k \ 2 + p^2wiW* k , (34) 

k=l k=l l^k 



with Lagrange multipliers X,p. Setting v := V] = VdFe 10 , equating to zero 

fc=i 

the derivative of (0) with respect to lOj and multiplying by Wj we get 

— \Wj\ 2 log I 1 2 + (A — l)|Wfc| 2 + p(v * Wj — \wj\ 2 ) = . 

Therefore, the quantity v*wjp and thus, after summing over j, also p, must be 
real, whence, necessarily Wj = e Vj, with Vj G R, for all j. The result follows 
by eliminating the overall phase. ■ 
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Proof (of Proposition 2.4): According to the previous Lemma, we choose \w) 
real and proceed to minimize 

d d d 

S(w,X,fi) := -J2 w l^wl + Xj2 w l + ^J2 w k ■ ( 35 ) 

k=l k=l k=l 

Because of convexity, the function g(x) := — x\ogx 2 intersects the straight line 
f(x) := 2(1 — X)x — /z in at most three points on [—1,1]. Therefore, the d 
solutions to 

— 2wk logw^ — 2wk + 2XiUk + /j, = , 

can have at most three different real values, a*, % = 1,2, 3. We denote by n» the 
number of times they appear among the components and consider the functional 

3 3 3 

S(a;nX,fi,u) : = - riia* log a* + A^r^a- + fi^Uidi , (36) 

i=l i=l i=l 

where we treat the 77^2 S clS continuous variables constrained by n\ + n 2 + n 3 = d. 
Minimizing ( |36"D yields the following equations 

rii(ai log a\ + ai — Xdi — /x) = , i — 1,2,3 (37) 
— log + Xa\ + /iOj + z/ , i — 1, 2, 3 . (38) 

It follows that, if n< > 0, i = 1, 2, 3, then, £- =1 (/zaj + 2z/ + 2a^) = 0, i = 1, 2, 3, 
and thus a = b = c. This case corresponds to Pf=i = I^X^I; a pure state, 
with null entanglement with respect to A. Therefore, if there are three different 
intersections, the minimum entanglement is reached at the boundary values of 
rij, 2 = 1,2, 3, that is, without loss of generality, at n\ — n 2 — 1 and n 3 = d — 2. 
If there are two intersections, that is if, without loss of generality, n 3 = and 
Qi 7^ 0,2 = «3, then, from (|37|j38D , we calculate /i = — 2(ai + 02), /i = 0102 and 
deduce the equality 



a 



a x — a 2 + a x a 2 log -2=0. 



a 



For fixed ax, because of their convexity properties, the two functions f(x) : = 

log — ^ and g>(x) := intersect at x = a 1; but, at no other points. Therefore, 

or a; a% 

the entanglement is again minimal at the boundary, that is at , say ri\ — 1 and 
n2 = d — 1. ■ 



Remark 2.4 Lagrange multipliers have been used in f|TT[| in order to calculate 



the entanglement of formation of isotropic states of bipartite quantum systems, 
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where it is shown that, when F > 1/d, the optimal decomposers have only 
two different components. We shall relate those results to ours in the following 



section, where we also discuss the fact, discovered in stating there is a 
bifurcation point F** such that the entanglement of formation is linear in F 
between F** and F = 1. 

Proposition 2.4 shows that when the vector \w) has only two different com- 
ponents, then we reduce to the case discussed in Proposition 2.3. Instead, when 
\w) has three different components, which is possible in a range of values of F, 
then we have more than one optimal decompositions. If d = 3 one gets at least 
two. Notice that these results are obtained under the hypothesis that Go-orbits 
of vectors \w) provide optimal decompositions for the entanglement of pp with 
respect to the subalgebra A. 

This fact is linked to the convexity of the function (|18|), which, as observed 
in the discussion of Case 2, fail in a neighborhood of F — 1: If F > F** one 
needs two orbits: the optimal orbit for F = F** and the singlet for F — 1, just 
as observed in fi~l|| . Consequently, for F** < F < 1 no Go-orbits can be optimal. 



3. ENTANGLEMENT AND ENTANGLEMENT OF FORMATION 

In this section we establish a one-to-one correspondence between the results 
of the previous section, in particular proposition 2.3, and the entanglement of 



formation of highly symmetric states as examined in [11 j. This concerns mainly 



the region (1/d) < F. From |T]J we learned the existence of the bifurcation 
point F**. On the other hand, our results in the region (1/d) < F < F** can 
be converted into those found by Terhal and Volbrecht. Indeed, the value of the 
entanglement of formation will be proved to be just fl3~0l) . 

To this end we consider the tensor product M. := A ® B of the full d x d 
matrix algebra, denoted by A, with a copy, B, of itself. We fix an orthonormal 
basis {\j)} of C d and given any density matrix, that is a state on A, 

PA = J2 R ^\j)(k\ , R=[R jk ]>0, TrR=l, (39) 

we embed it as -D^a] into the state space of A4 according to the following 
Definition 3.1 Let D be the linear map associating matrix units \j)(k\ of A 
with matrix units <8> \j)(k\ of M.. We shall refer to it as the doubling map. 

It transforms states pa on A into states on Ai = A <S> B of the form 



14 



PA ^ D[p A ] :=Y,Rjk\j)(k\ ® \j)(k\ , (40) 

Remark 3.1 This yields the class of density matrices in Example 2.3, which 
we shall refer to as diagonal class (with respect to the chosen basis). On the 
given diagonal class the doubling map can be inverted 

D- 1 : Pab = J2 R i,k\j)(k\ ® \j)(k\ ^ p A =Y, R iM(k\ ■ (41) 

j,k j,k 

The argument developed in Example 2.3 ensures that decompositions of pa can 
be mapped onto decompositions of -D[p^]. Vice versa, decompositions of pab 
provide decompositions for the diagonal class of pa by applying D~ x . Moreover, 
if Aq C A denotes the subalgebra of diagonal matrices in the given, fixed 
representation, then S(p \Aq) = S(D[pa] \A). Therefore: The entanglement is 
preserved by D, in the sense that 

E(p A ; A, Ao) = E{D[ PA ); A ® B, A) . (42) 



In [[Ll|] the entanglement of formation has been calculated for the isotropic 

states 

u f = ^^(1ab-\*)(*\)+F\*)(*\ . (43) 
In the above expression Iab is the identity for the algebra A <S> B and 

l*> = 4j£L7>®b'> • (44) 

V« j= i 

Remark 3.2 The isotropic states are invariant under the group Q of all uni- 
taries of the form U ®U where (a\U\b) = (a\U\b)*, 

U ® UujfU- 1 ® (J' 1 = cu F . (45) 

As in Remark 2.2, it follows that u F can be expressed as the following average 
with respect to the Haar measure dcU, 

u F = [ dglJU ^Ulty^lU- 1 , (46) 

Jg 

if and only if F = {9f\u F \^f) = |(*|$)| 2 . 

We compare the isotropic state u F with the doubling of p F in (PH), 
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D ^ = J^j( D lU-Dm(i>}\) + fd 



= 73t(E ® \m - l*X*]l) + *1*X*I • ( 4 7) 

« 1 j=l 

Proposition 3.1 Let F > 1/d and consider the decomposition 

by means of the unitaries £/„- that implement the permutation group G. If the 
latter is optimal for the entanglement of formation E(up) with |$)($| in the 
diagonal space, then E{uj f ) = E(p F , A, Ao). 

Proof: The d\ unitaries form a subgroup G <S> G of the group of unitaries 
in Remark 3.2; they implement the permutation of the chosen basis {\j) ® \ j)} 
of the diagonal space. Then, (\I/|u; j p|\]/) = (^|£)[pir]|^) = F and 

d[pf] = I E^ 1 ® tf,rW ® £4 • 

If |$)($| is optimal for up, it turns out from Proposition 2.2 that the decom- 
poseres U <8> t^l^X^t/ -1 ® t/ -1 are optimal, too. Thus the result follows from 
Proposition 2.1. ■ 

Remarks 3.3 

(i) If F > 1/d the isotropic state u F is entangled. When F < 1/d it becomes 



separable. There exist several proofs of this fact, e.g. [18 



(ii) In view of Remark 2.3(h), the previous proposition establishes a link be- 
tween our results and those of |LlJ| . In |ll] a new symmetry breaking bifurcation 



point was observed at F = 8/9 when d = 3. The doubling map makes it corre- 
spond to a bifurcation point within case 2 of the previous section at the same 



value of F The numerical analysis in |T4| missed it, the needed accuracy being 
of the order of 10~ 4 . In both cases the leaves TZ(u F ), respectively TZ(p F ), are 
identical for all F within F** = 8/9 < F < 1. This unique leaf is generated by 
the optimal decompositions of u 8 / 9 respectively p 8 / 9 , which form one orbit, and 
by the pure state u\ given by ( fPj ) respectively p\. The latter orbits are singlets. 
(Hi) The entanglement of p\ and p 8 / 9 that generate the leaf discussed in the 
previous remark do not coincide, 

1 

3 



E( Pl ;M,A) =ln3 , E(p 8/Q ;M,A) = In 3 - -In 2 . (48) 
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We shall now relate the remark above to another observation which again 
relate entanglement of different algebras with one another. 



From Case 1 in section 2, we know that vectors of the form 



and 



y . 



y 

X 



with x 2 + y 2 = 1 generate the leaf of some state p 2 on M 2 (C). These 



2-dimensional vectors can be embedded in C 3 as follows, 





( X \ 




( y \ 


\wi) = 


y/V2 


, k 2 ) = 


x/V2 




\y/V2j 







(49) 



With them we construct the density matrix in M 3 (C) of the form 
p 3 = X\wi){w 1 \ + (1 - X)\w 2 ){w 2 \ 



fa 


b 


b\ 


b 


c 


c 


u 


c 





(50) 



It is easy to check that powers of p 3 have the same structure which is thus in- 

fu\ 

herited by yfps . It thus follows that \fp%\(f)) = 



v 



for any The discussion 



of Example 2.3 assures and that the optimal decomposers of p% with respect to 
the entanglement E(p; M^(C) , A3) , with A3 the maximally Abelian subalgebra 

in the chosen representation, have again the same form. But then, being 



x 



and 



y 

X 



optimal with respect to E(p 2 ; M 2 (C), A2), QSTID is itself an optimal 

decomposition of p 3 with respect to E(p 3 ; M 3 (C), A3). 

According to the discussion at the beginning of this section, it also follows 
that the doubling map 



\wi) 
\w 2 ) 



\Wi) = x\l) ® |1) + ^(l 2 ) ® |2> + |3) ® |3 



W 2 ) = y\l) ® |1) + -j=(\2) ® |2) + |3> ® |3)) 



(51) 
(52) 



provides optimal decomposers, too. In particular, for given x, y on the unit 
circle the pure states j = 1, 2, generate a leaf of the entanglement of 

formation functional on which it is convexly linear. 

Moreover, for x = l/y/3 and y = y2/3, we get \Wx) = |^), with fidelity 
F = K^l^i)! 2 = 1, and \W 2 ) = |$ 8/9 ) with fidelity F = \(V\W 2 )\ 2 = 8/9, 
indicating a reason for the bifurcation value F = 8/9. 
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One observes that (51) and (|52[) become identical for x — y — 1/ v^2 so that 
the doubling map gets the vector 

\Ws) = i|l> ® |1) + i(|2> ® |2) + |3) ® |3>) , (53) 

which has fidelity 

FH(*l^3>| 2 = ~ + yf = P+(l-p)~ , 0<p = 3v / 6-^<l . (54) 

Let us now consider the state 

p F = j?|tf>(*| + (1 -p)|$ 8/9 )($ 8/9 | . (55) 

By using (|4"8"D, it can be shown that its entanglement E(pp) is larger than 
pE(p{l)) + (1 -p)E(p(&/9)) for < p < 1. This implies that convexity of 5(F) 
in (|3~2"D is lost for F > F** in accordance with the discussion above. 

We finally note that one can extend ( J4"9| ) to all dimensions larger than two. 
Indeed, let z±, z 2 denote the components of a unit vector in two dimensions. By 
similar arguments one proves that the leaves of case 1 of the previous section 
are mapped onto certain leaves belonging to the entanglement of formation in 
d + 1 dimensions by the embeddings 

( Zl ) — zi\00) + (z 2 /Vd) 4 f:\jj) (56) 

In particular, the embeddings of {zi, z 2 } and {z%, zl} form an optimal pair with 
respect to the entanglement of formation. One further observes in the special 
case Zi = 1/ Vfl+1 the embeddings (|56|) are the totally symmetric vector \I/ in 
d + 1 dimensions and 



d , , / 1 ^, , , , 



V tZ+l V l ) j: 1 

Its fidelity reads F = 4d/(l + d) 2 , and we see as above 

Fr +1 =Hd+l)- 2 (5 



i. e. the bifurcation value given in [IT] for d + 1 > 2. 



4. CONCLUSIONS 

We have studied in several examples the entanglement defined by a maximal 
commuting subalgebra of a full matrix algebra, and in its relation to the entan- 
glement of formation. Apart from its actual numerical value, what is interesting 
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is the structure of both entanglement functionals upon the space of states, and 
their separation into different leaves. To some extent these leaves can be found 
by applying group theoretical considerations. They show a rich structure with 
varying stability under the groups under consideration, Since the same group 
appears in different algebraic contexts, it can be shown that the decompositions 
of states on different algebras can be related. This helps to control the optimal 
decompositions and to understand their variety. This new technique is shown 
at work in several examples: The doubling map relates two quite different lines 
of research which had been considered almost independently up to now. In par- 
ticular we have a further proof of the entanglement of formation results for 



isotropic states of Terhal and Volbrecht in the region (1/n) < F < F**, [11 



Another embedding map verifies their bifurcation point F** close to F = 1 as 
a footprint of a symmetry-breaking in two dimensions. It belongs to class of 
maps which change entanglement but not the leaves. The leaves are respected 
because the entanglements differ just by a convexly linear function. 

It should be clear that we only provide some distinguished first examples of 
our embedding procedures which can connect various entanglement problems 
and, evidently, other ones which are defined via convex or concave roofs, for 
example general entanglement monotones or Holevo (1-shot) capacities. 
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